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1 Introduction 



This paper results from a study of a control problem for a thermodynamic process gov- 
erned by the reaction-diffusion equation: 



Ut{x, t) — A'u(x, t) = f{u{x, t)) 



The underlying idea refers to an implementation of the closed-loop control system with 
a finite number of both control units and measurement points inside the domain of the 
process evolution. Such a control system leads to the following system of parabolic partial 
differential equation coupled with ordinary differential inclusions: 



Ut{x, t) — Au{x, t) = f{u{x, t)) + g{x, t; K(t)) 
f3iki{t) + Ki{t) e Wi{t,ulxl,t),...,u{x*„,t)) 



on 
on 



Qt 
[0,T] 



, Pmfimit) + f^m{t) ^Wm{t,u{xl,t),...,u{xl,t)) OU [0, T] 



;i.2) 



du 
dn 



u{0, x) = Uq{x) 



on on X (0, T) 
for X G 

for j = 1, . . . ,m 



where Qj- = Q x [0,T), Q is bounded, open and sufficiently regular domain in M'^, x^. G Q, 
k = 1, . . . ,n are fixed points in the interior of Q, Wj defined on [0, T] x MJ^ for j = 1, . . . , m 
are given multivalued functions, / : M — )■ M and g : x — > M. The unknown is 
represented by the couple {u, k), where u : Qt — i^j '■ [0,T] — )■ M and K,{t) = {K,j{t))^^^. 
The main objective of the present paper is to explore existence questions for the 



solutions to system (1.2). 



The idea of controlling a thermodynamic process by a closed-loop control based on a 
finite number of control devices was already addressed in several publications. To mention 
some of those, problem of controlling linear heat flow was considered for example in [7j 
(for one-dimensional model) and in j8], [9], [lO] (for multi-dimensional models). Closed- 
loop controls of the above type for the two-phase Stefan problem was also extensively 
investigated, cf. [6j (one-dimensional model) and [llj . ^ (multi-dimensional case). A 
nonlinear parabolic system without free boundary was treated in j3] in control context, 
however a very special form of the nonlinear term was admitted there. 

In general, the concepts of closed-loop control systems presented in the above papers 
can be divided into two categories: 1) the control and measurement devices are supposed 
to take their actions without delay and 2) the occurrence of such delay is admitted. 



The model described by (1.2) corresponds to the former type. We also try to keep 



assumptions concerning the nonlinear term / possibly general. Moreover, the structure 



of (1.2) assumes that the control devices are placed inside the domain of the process (the 
term g represents the control devices actions), while in all above mentioned papers the 
control devices were assumed to act through the boundary of the domain. 

The plan of this paper is as follows. In Section [2] solutions to system (1.2) are defined 
and an existence theorem for those solutions is formulated. 
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Section [3j shows an applied context for the system (1.2): a construction of the closed- 
loop control system based on a finite number of measurement points and a finite number 
of control devices is presented. 

Section |4] provides a proof of the existence theorem stated in Section |2} The proof is 
based on a fixed point argument for a multivalued operator with the use of a generalized 
Kakutani fixed point theorem. 



2 Definition of solutions for system (1.2) and the main 
existence theorem 

We first formulate a set of basic notations and assumptions for this paper. 
Let us denote 

K = C{[0,T];W) (2.1) 

and, for any arbitrary R> 0, 



M, 



R= {(«i,...,M G W^''°°(0,T;M'") : || ||^i,^(o_^) < i?, z = l,...,m} . (2.2) 

Mji is nonempty, convex and compact subset of with compactness following via the 
Arzela-Ascoli theorem. 
We assume that: 

(I) VL is an open and bounded domain in M*^, with C^-boundary. 
(II) G int for A; = 1, . . . , n. 

(III) a) / satisfies the linear growth condition, f{s)s < Ci +C2S^ for some positive finite 

constants Ci, C2, 

b) / is locally Lipschitz continuous, 

c) /(O) = 0. 

(IV) The multivalued functions Wj : [0, T] x — > 2^, j = 1, . . . ,m satisfy the condi- 
tions: 

a) Wj{t,r) are nonempty, closed and convex subsets of M for every (t,r) G [0, T] x 



b) Wj are upper semicontinuous (see Definition 5.1), 

c) Wj are bounded: there exist finite Cj > such that W{t,r) C [—Cj,Cj] for 
every (t, r) G [0,T] x M". 

(V) We set: 

a) S = maxj=i,...,m |l «i I + Cj + ' ""^'"^^ | , 
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b) p and q are fixed reals such that 

p > 1 and q >2 and 1 — < 1 

2p q 

in particular, q = 2 and p = d. 

(VI) a) ^(., .;/?(.)) e L2(0,T;L°°(fi)) for every K G ir, 

b) \\g{., - g{., llLi(Qr) ^ CsEJli Ik' - «^ 



'i i llc[o,T] 



with some positive constant C3, for all k ,k G Ms 

c) II fi'l • ) • ; • )) II L'J(o,r;LP(n)) — '^4 

with some positive constant C4, for every k G Mg. 
(VII) Mo e L^{n). 
(VIII) /3j > 0, G M. 



Prior to defining solutions to the system (1.2), we shall define those of the following 
parabolic problem: 

Ut{x,t) — Au{x,t) = f{u{x,t)) + g{x,t) on Qt 

g = on dnx (0,T) (2.3) 

u{0,x)=Uo{x) for X E fl 



and of the ordinary differential inclusion 



/37(t) + lit) G W(t) on [0, T] 
7(0) = 70 



(2.4) 



the latter representing a generalization of 



/37(t) + 7(t) = V(t) on [0,r] 
7(0) = 7o 

For the system ( |2.3[ ) we assume the following definition: 

Definition 2.1 [] For Uo G L^(n) and (7 G L'^{0,T; H^(Q,)*), u is a weak solution to the 
problem (2.3) if 

u G L2(0,T;ifi(l])) nC([0,T];L2(fi)), v! G ^^(0, T; ifi(l])*) 

and ^ 

[ {(j),u')~{(f)J{u)+g) + {V(j),Vu)dt = 
Jo 

for all (p G L2(0, T; ^^^(^^)) and u(0) = no. 



^(., .) denotes the paring between H^{n) and its dual space, (., .) denotes the standard scalar 
product in L^{yi). 
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For g is of the form g{. , . ; /t( . )), meeting the assumption (VI) and k, ^ K, the problem 
(2.3) admits a unique weak solution u (see Theorem 5.1) which is continuous in Qx (see 
Theorem 5.3), thus a (single- valued) mapping 

u:K^ C{Qt) (2.6) 



can be defined, which to any given n G K assigns the weak solution u of (2.3). 
The solutions of (2.4) are defined by the following 

Definition 2.2 For 70 G M an(^ W G L°°(0, T): 



a) "J is a solution to (2.4), if in the Caratheodory sense it is a solution to (2.5) for some 
integrable function V : [0,T] — >R such that Y{t) G W(t) for a.e. t G [0,T]. 



b) is a solution to (2.5) in the Caratheodory sense, 2/7 G AC[0,T], 7(0) = 70 and the 
first equation in (2.5) is fulfilled by 7 a.e. on [0,T]. 



Remark. The above definition remains correct also for a wider class of functions 
W, still since in assumption (IVc) we put a boundedness constraint for Wj appearing in 
(1.2) it appears convenient to focus on the case of W G L°°(0,T). 

If given functions Wj, j = 1, . . . ,m satisfy the assumption (IV) and u G C{Qt), then 
for Wj = Wj{t, u{xl,t), . . . , m(x* , t)) entering into the right-hand side of (2.4) there exists 
a solution in the sense of Definition 2.2 (see Theorem 5.5) what allows us to define, for 

(2.7) 



every Wj, j = 1, . . . , m, a multivalued mapping 



which assigns the set of solutions of (2.4) to given u G C{Qt). 

Proposition 2.1 LetW be given such that a solution to the equation ( 2.4^ (in the sense 
of Definition 2.2) exists. Then 'y is a solution to (2.4) if and only if 

7(t) = e~?*7(0) + ^ / e--^^'-'^Y{s)ds 
P Jo 

for some integrable function V : [0, T] — > R such that V(t) G W(t) for a.e. t G [0, T]. 

We are now ready for formulating the following definition: 

Definition 2.3 Suppose that il, f, g, Wj, Uq, aj and 13 j satisfy assumptions - (VIII). 
We say that a couple {u,k,) G C{Qt) x [AC[0,T]]'" is a solution to system (1.2) if and 
only if 



a) u = u(ft), i.e. u is a weak solution of system (2.5) in the sense of Definition 2.1 
corresponding to k. 
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b) Kj G ij{u) for j = 1, . . . ,m, i.e. is a solution in the sense of Definition 2.2 of (2.4) 
with Wj{t,u{xl,t), . . . ,u{x*,t)) entering into the right-hand side, or equivalently (by 



Proposition \2. 1 ) 



K,{t)=e-^^\M + ^ f e-^^'-'\,{s)ds 

Pj Jo 

for some measurable function Vj : [0,T] — )■ R, such that 

vj{t) e Wj{t,u{xl,t),...,u{x;,t)) 
for a.e. t G [0, T] for j = 1, . . . ,m. 

We conclude this section with the main result of our paper: 



Theorem 2.1 Suppose that assumptions as stated in Definition \2. 3 are fulfilled. Then 



there exists a solution to system (1.2) in sense of this definition. 
The proof is given in Section |4] 

3 Construction of the closed-loop control system 

As mentioned in Introduction, we assume that the control of the diffusion-reaction process 



( 1.2 ) is performed via a finite number of control devices localized in the interior of the pro- 
cess domain. The control action is based on observation data at finite set of measurement 
units as compared to given reference patterns. 



Suppose that a reaction-diffusion process governed by (1.1) with homogeneous Neu- 
mann boundary conditions and prescribed initial condition evolves over Q. The control 
action for this process is then described by an additive control term (denote it by g) to 



the reaction-diffusion equation assuming so the form (2.3). 

Let the control term be set as g = g{x, t,C,), where ^ G is an m-dimensional param- 
eter. Moreover if, for any admissible function k, : [0,T] — > M™, the term g{x,t, K{t)) is 
sufficiently regular, then any such n can be considered as a control applied to the system 



(2.3) and it determines a unique solution u which represents the response of the system 



to K. 

To be more specific, it is our aim to keep the state of the system „as close as possible" to 
a prescribed function it : fl — )■ M which assumes values ul := u{xl) at the measurement 
points. 

Let us expose more on the above idea of measurement points and control devices. The 
control K, shall be synthesized by m control devices located in int Q which have an ability 
to drive the thermodynamic process in the entire domain: each of m components of k 
is considered to be an input signal for one of those m devices. We implement the above 



concept by assuming the control term g in (2.3) to admit the representation 



g{x,t,K{t)) = ^gj{x,t)Kj{t) for {x,t) G Qt , (3.1) 
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where gj are prescribed nonnegative functions on Qt- Each of the functions gj describes a 
[ki, . . . , Km) and Kj is the signal which specifies actions of j-th device. 



control device, k 

described as the solutions of 



K,j{0) = Qj 



on 



[0,T] 



(3.2) 



for j 



[0,T] 



m, where Pj are positive constants characterizing the j-th device and 
I can be considered as an external power supply term for j-th device. 



Now let us assume that for controlling the behaviour of system (2.3) we use a control 
K which itself depends on a solution u of (2.3). To be more specific, we postulate that the 
latter dependence is included in the power supplies Vj and admits the representation 



for j = 1, 



W = X] (^jk{t)wk{u{xl, t) - ul) on [0, T] 

k=l 

1, . . . , n, where ajk satisfies: 

a,k e C[0, T] 
a,k{t)>0 on[0,r] 
ELi '^jkit) = 1 on [0, T] for j 



(3.3) 



(3.4) 



1, 



, m 



Here Vj are convex combinations of Wk{u{xl, t)—ul.)) with t-dependent weights. It remains 
to define the functions w^- Through many concepts concerning Wk, the simplest possible 
idea, also of practical relevance, is alternating their values between —1 and +1, depending 
on the sign of u — ul at the control units xp 



Wk(r 



-1 


for 


r 


> 





for 


r 


= 


+1 


for 


r 


< 



(3.5) 



Taking into account ( 3.1[ ) - (3.5), we obtain a new system of differential equations with 
the unknown the couple {u, k): 



' ut{x,t) - Au{x,t) = f{u{x,t)) + Yl%igj{,x,i)Kj{t) 
Piki{t) + Ki(t) = Yl=i(^ik{i)wk{u{xl,t) - ul) 

I3m^m{t) + Km{t) = YJk=l Olmk{t)Wk{u{xl,t) - ul) 



on Qt 
on [0,T] 



(3.6) 



on 



[0,T] 



subject to suitable boundary and initial conditions. This is an algebraic system, in the 
form resembling (1.2). Still, the proof of Theorem 2.1 presented in Section |4] shows that 
assumption (IVb) on the upper semicontinuity of Wj entering into (1.2) is essential. Hence 
our existence result does not apply directly to (3.6) where Ylli=i^jkif)'^k{u{xl^t) — ul) 
stands for Wj and Wk are as in (3.5). This suggests an idea of replacing Wk with their 
multivalued convexifications Wk which are upper semicontinuous: 



Wk[r) 



+1 



for r > 
for r = 
for r < 
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and, accordingly, replacing the ordinary differential equations in (3.6) by the differential 
inclusions: 



Pjkj{t) + K.j{t) e ^ajk{t)wk{u{xl,t) - ul) 



k=l 



on [0, T] for j 



m. 



Now note that assumption (IV) is fulfilled for Wj = Ymc=i '^jk{t)wk{u{xl.,t) — ul) and. 



after ensuring that the rest of assumptions (jlj) - (VIII) is satisfied, including the choice of 
Qj such that g satisfies assumption (VI), Theorem 2.1 becomes applicable. This justifies 



the use of the differential inclusions in (1.2), since it seems natural to expect that the case 
of the switching control law inscribed in Wj will be taken into account. 



4 The proof of the existence theorem 



Let us now proceed to the proof of Theorem 2.1 The proof will exploit a fixed point 



argument. More precisely, it will be shown that if a certain multivalued operator on a 



Banach space has a fixed point (see Definition 5.3), then this fixed point determines a 



solution of problem (1.2) in the sense of Definition 2.3 (see Subsection 4.1). Next we will 
show that the discussed operator indeed admits at least one fixed point (see Subsections 



4.2 and 4.3) 



4.1 Construction of the fixed-point operator 

This subsection is devoted to constructing the suitable multivalued operator whose fixed 



points determine the solutions of problem (|1.2|) in the sense of Definition 2.3 

= 1, 



Let K be defined as in (2.1), Mr be defined as in (2.2), x 



1, . . . , n, meet 

assumption (IIT]) and S fit assumption (Va). For every ^ G M5, the solution u = u(^) is 



continuous on Qt (where u is the operator defined in (2.6)), hence the functions u{x*f,, 
1, . . . ,n are well-defined continuous functions on (0, T). Therefore the operator 



k 



R:Ms^C{[0,T];. 
R(0 = (u(0(xt, .),. 



.u(e)«, .)) 



is well-defined. Referring to the ideas introduced in Section |3j operator R incorporates 
the whole information about the process at given control points. 

Denote W = {Wi, . . . , Wm). We define a multivalued operator Q : C([0, T]; M") — > 

for given / G C([0, T]; M"), any v G L°°(0, T; R"") satisfies v G Q(/) 
if and only if v{t) G W{t, f{t)) a.e. on (0,T) . 



By assumption (IVc), each Wj is bounded with some finite constant Cj > 0. Every 
V = (fi, ...,Vm) e L°°(0,r;M™) that belongs to U^^(Q) satisfies 
j = 1, . . . , m. Denote C = {Ci, . . . , Cm) and 



"^j IIl°°(o,t) — 



c 



V 



[Vl, 



.Vr. 



G L°°(0,T; 



V 



j IIl°°(o,t) 



< C, for J = 1, 



, m 
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The set Mc is bounded in L°°(0, T; M™). As in Section |3| the operator Q determines a 
control law which assigns the power to be supphed for the control devices to information 
about the process collected in measurement points k = 1, . . . ,n. 

The values Vj determine responses of the control devices kj, j = 1, . . . ,m described 
by ordinary differential equations 

[],k,it) + K,it) = v,it) on [0,T] 

the latter having solutions (in the Caratheodory sense) given by the integral formula 
introduced in Proposition 2.1| with Vj instead of V: 



Kj(t) = e a,- + T^e-*/^^ / e'/^^ Vj{s) ds . (4.1) 
Pj Jo 

This formula sets continuous functions well-defined on [0,T] for any Vj G L°°(0, T; M"^): 



thus we can define an operator P : L°°(0, T; M"^) — > K with {P{v))- given by (4.1) for 
j = 1, . . . ,m. As in Section |3j the operator P assigns an input signal for control devices 
corresponding to the supplied amounts of power. 

Moreover, the solutions kj are absolutely continuous and bounded, with bounded 
derivatives, hence kj G W^'°°{0,T). We can also conclude the following estimates for the 
values of kj and kj under assumption that v G Mc'- 



Jo I3j ' 







= I a, I + (1 - e-*/^0 < |a,|+C, 

for j = 1, . . . , m. In the same way we show that Hj{t) > — \aj \ — Cj. For estimating kj, 
j = 1, . . . , m we proceed as follows: 

—Cj < /3jkj + Kj = Vj < Cj , 

thus 



Ki ^ ^ < ^ ^ 

since K,j{t) G [— | o-j | — Cj, \aj\ + Cj]. Similarly we obtain kj > _\jhhL^_ Altogether for 
all J = 1 , . . . , m we have 

aj\ + 2 



II vKi.°°(o,T) — j^^'^q 1 1 '^i I ' ^ 



Corollary 4.1 For every v G Mc, P{v) G Ms where S is as in the assumption (Va) 
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Actually, the above corollary justifies the choice of the constant S introduced in assump- 
tion (fVal). 



To conclude, we have defined the following operators: 



R 

Q 
P 



M5^C([0,T];M") 



C([0,T];: 
L°°(0,T;: 



2-L°°{0,T;M™') 

K 



(4.2) 



such that 



: Mr 



\Mc ■ ^'^c — ^ Ms, 

thus the superposition P o Q o R : M5 — > 2^^^ is a well-defined multivalued operator. 

Note that by definitions of the operators P, Q and R, the condition K G P o Q o R(e) 
can be rewritten to the form: 



1 



e Vj{s) ds for j 



m 



for some measurable function Vj : [0,T] — y R such that Vj{t) G Vrj(t,R(^)) for a.e. 
t G [0,r], where iy,(t,R(0) = t), . . . ,M(x;,t)) on [0,T] for m = u(0- Equiva- 

lently, in more compact form: kj G tj{u) for u = u{^) and j = 1, . . . , m (where tj are the 
operators defined in (2.7)). 

This brings us to the conclusion that for ^ = k the condition k G P o Q o R(^) is 
equivalent to the hypothesis included in Definition 2.3 Thus, we can reformulate the 
definition of solutions to system (1.2) by introducting the following equivalent: 

Definition 4.1 Suppose that Q, f , g, Wj, Uq, aj and (5j are obeying general assumptions 
say that a couple {u, k) G C{Qt) x 74C([0, T]; M™) is called a solution to 
system (1.2) if and only if 



a) u = u(/t), i.e. u is a weak solution of system (2.5) in the sense of Definition 2.1 
corresponding to n, 



b) ne PoQoR(k) 



The latter definition allows us to conclude that if k is a fixed point of the operator PoQoR, 
i.e. K G P o Q o R(/t), then the couple (u(k), k) is a solution of (1.2), hence the existence 
of solutions (1.2) will be assured once we prove the existence of at least one fixed point 
of P o Q o R. The existence of the fixed points will be shown in the next two subsections. 

Remark. By the construction of operator P o Q o R we conclude that solutions 



{u,k) to system ^L2^ belong not only to C{Qt) x AC{[0,T]; 
PF^'~(0,T;M'"). 



but also to C{Qt) 
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4.2 The fixed-point argument 

As already mentioned, we shall prove the existence of fixed points for the operator PoQoR 



defined in Subsection 4.1 To this purpose, we will employ a certain generalization of the 



Kakutani fixed point theorem — we will take three abstract operators Q and IZ (not 



to be confused with P, Q and R) acting between spaces given in (4.2) and investigate 



their properties that facilitate an appropriate use of the generalized Kakutani theorem. 



Once we give an answer to the above problem, we pass in Subsection 4.3 to check- 



ing that operators P, Q and R defined in Subsection 4.1 indeed have got the required 
properties established for V, Q and IZ. 

To start our considerations, we recall the following generalized Kakutani fixed point 
theorem ([21 Theorem 4]): 

Theorem 4.1 Suppose that X is a real Banach space and let M G X be nonempty, 
convex and compact. Suppose also that the multivalued operator T : M — > 2^^ satisfies 
the conditions: 

a) T has nonempty, closed and convex values (i.e Tx has these properties for all x & M), 

b) T has closed graph in X x X . 

Then there exists an element x E M such that x G Tx. 

□ 

Before proceeding further, for a sequence (/*)^i of functions in C([0, T]; M") (denote 
fi = {fl . . . , /;)) and a function / G C([0, T]; M") (denote / = (/i, . . . , /„)), we introduce 
following property which will be useful in the sequel: 



for all t G (0, T), S > there exists io{t, 5) G N, ho{S) G (0, S] such that 
max \ fl(t + h)-fk{t)\<6 (4.3) 

for all i > io{t, 6),0<\h\< ho{S),t + h e {0,T) . 



Now, let K be defined as in (2.1) and be defined as in (2.2). Suppose also that 
a multivalued operator Q : Mr is of the form Q = V o Q o TZ where V, Q and V 

satisfy: 



Hypothesis 4.1 The operator TZ satisfies: 

a) n:Mn^C{[0,T];R^), 

b) for any sequence k in Mn with topology of K, the sequence /* = TZ{k,^) has a 



subsequence convergent to f in sense of (4-3) with f = TZ{k). 



Hypothesis 4.2 The multivalued operator Q satisfies: 
a) Q : C([0,T];M") -^2^°°(0'^'K'"), 
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b) Q is closed in the following sense: if 



'>') (Dill C([0,r];]R") converges to f in sense of {4-3), 
ii) v' m L°^(0,T;M"^), 
Hi) G 

then V e Q{f), 

c) Q is bounded, i.e. [JIm{Q) is a bounded subset 0/ L°°(0, T; M™), where Im{Q) :- 
Q(C([0,T];M")). 



Hypothesis 4.3 The operator V satisfies: 

a) V : L°°(0,T;M'") — > K, 

b) ^lu/m(S) : [jIm{Q)-^MR, 

^) im{Q) is closed in the following sense: if 

i) V in |J/m(Q) with topology of L°°{0,T; 

ii) K,"^ ^ K in Mji with topology of K, 

lii) = P|u/m(Q)(f*) 

then K = V\\jim{Q){v). 



The following lemma holds: 



Lemma 4.1 Under Hypothesis \4.l\ \4-^ and\4_^ Q = V o QoTZ : Mr 

operator, with topology of K used in Mr. 



2Mr 

is a closed 



Proof. Suppose that ^ k and ^ ^'m Mr with topology of K and that G Q{k^). 
The question remains open whether ^ G Q{k). 

It is equivalent to say that ^ = V{v) and v E Qo TZ{k) for some v G L°°(0, T; M™). 
At the same time, ^* G Q{k'^) is equivalent to ^* = and G Q o TZ{k'^) for some 

V G L°°(0, T; M™). Due to Hypothesis 4.2, is a bounded sequence, thus there exists 



a weakly-star convergent subsequence (for simplicity, we keep the original indexes when 
passing to a subsequence): 

A ^ in L°°(0,T;M'") . 



By Hypothesis 4.3, we hence obtain = V{v) 



Now it remains to check that f G Q o TI{k), i.e. v G Q{f) and / G TI{k) for some 
/ G C([0, T]; M*"). At the same time, the condition 1;* G Q o TZ[k^) is equivalent to 

G Q(f ) and f G 7^(/^0 for some f G C([0,T];M"). Due to Hypothesis 
exists an element / and a subsequence of /* such that (again after reindexing) 



4.1 



there 



/*-)■/ in sense of (4.3) 



and moreover / = TZ{k,). By the above property, weak-star convergence of and Hy- 



pothesis 4.2, we can claim that v = Q{f)- 
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To sum up, we infer that there are elements / and v such that / = 'R-{k), v = Q{f) 
and ^ = V{v), the latter completing the proof. 



□ 



The next lemma follows straight forward as the corollary of the former one: 



Lemma 4.2 Under Hypothesis 4-1^ 4-^ and 4-3, the operatorQ = VoQoTl : Mji — > 2^^ 
has closed values. 

□ 



Aiming to fulfill the assumptions of Theorem 4J^, we still need to guarantee that Q 
has nonempty and convex values. We need to make another hypotheses which will assure 
these properties: 

Hypothesis 4.4 V is an affine operator, i.e. for every x G L°°(0, T; M™) we have V{x) = 
Ho + "Poix) where yo ^ K and Vq : L°°{0, T; M."^) — > K is a linear operator. 

Hypothesis 4.5 Suppose that a multivalued function W : [0,T] x M" — y 2"^™, W = 
{Wi, . . . ,Wm) is given where Wj satisfy the general assumption (IV). We assume that Q 
is defined in the following way: 

for given f G C([0,T];M") any v G L°°(0,T;M™) satisfies v G Q(/) if and only if v{t) G 
W{tJ{t)) a.e. on (0,T). 

Lemma 4.3 Under Hypothesis \4-^ the operator Q: 

a) has nonempty and convex values 

b) is bounded, i.e. there is come C > such that Q{f) C i?ioo(o,T_R™)(0, C) for every 
/gC([0,T];M") 

Proof. By the boundedness of W, the second assertion follows directly. The convexity 
also is in turn a straightforward consequence of the convexity of the values of W. 

It remains to prove that for arbitrary / G C([0, T]; M") there is some v = [vi, . . . , Vm) ^ 
L°°(0, T; M™) which belongs to or, equivalently that Vj{t) G Wj{t, f{t)) a.e. on (0, T) 

for j = 1, . . . ,m. By the boundedness of W, it is enough to check that a measurable 
selection is possible, namely that there exists some measurable vj satisfying the last 



condition. But this is assured by Theorem 5.4 



□ 

It follows from the above result that: 

Lemma 4.4 Under Hypothesis \4-i\ U-^j \4-4\ and \4.5[ the operator 
Q = V o QoTZ : Mfi — y 2^^^ has nonempty and convex values. 
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Proof. The non-emptiness is a consequence of Lemma 4^ and the well-posedness of 
the operators V and TZ. 

Let us still check the convexity. Suppose that G Q{k) for some k G Mr. We 

want to show that for arbitrary A G (0, 1) there holds := A^^ + (1 — A)^^ G Q{k). It is 
equivalent to say that = V{v^) and f G Q o 71[k) for some v^. 

Note that the condition C,^ G is equivalent to = V{v^) and f * G Q o 7^(/t) for 



some V 



1, 2. Thus let = + (1 - XW- 



v"- G 



Q o 7^(fi;) because Q has convex 
values, as stated in Lemma 4.3 Now, taking into account the affinity of ^ = V{v^) 



can be obtained by direct calculation. 



□ 



Due to the above lemmas, under Hypothesis 4.1 - 4.5 it is possible to use Theorem 4.1 



still there is one more one simplification that can be made, namely it can be shown that 
Hypothesis 4.5| is stronger that Hypothesis 4.2 



Lemma 4.5 If the operator Q satisfies Hypothesis \4-5\ then it also satisfies Hypothesis 
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Proof. Condition afi in Hypothesis 4.2 follows directrly, condition ch is assured by 



Lemma 4.3 The proof that the condition |b]) is assured can be found in P, theorem 3.6] 
for the case n = m, an analogue of this assertion was used in [TT| lemma 4.4] for n ^ m. 

□ 

Altogether, as a concluding corollary we can state the final theorem of the above consid- 
erations: 



Theorem 4.2 Suppose that Q = V o QoTZ where V , Q, TZ satisfy Hypothesis 4-i\ \4-^ 
4-4 and \4.5[ respectively. Then there exists a fixed point of Q in Mr, i.e. there exists a 
K G Mr such that k, G 



Proof. We can apply Theorem 



assured by Lemma 4.5, Lemma 4.1 



4.1| with X = K and M = Mr. Its hypotheses are 



(closedness of the fixed-point operator). Lemma 4.2 



(closedness of its values) and Lemma |4.4 



□ 



Remark. The proofs of Lemma and AA exploited some arguments utilized in the 
proof of [11, theorem 4.1]. 



4.3 Conclusion of the proof of existence 

This subsection is devoted to checking that the operators P, Q and R from Subsection 
4.1 1 satisfy Hypothesis 4J^, 4^, 4A and 4^ stated in Subsection 4^ As stated in Theorem 
4.2| this will imply that the superposition P o Q o R has at least one fixed point k, G Ms 



(where Mr is defined in (2.2) and S is the constant set in assumption (Va)). Due to the 
Definition 4.1 (equivalent to Definition 2.3), it will mean that the fixed point k, determines 
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a solution {u{k,), k) of the problem (1.2) (where u is the operator defined in (2.6)), hence 



completing the proof of Theorem |2.1[ We remind that for Theorem 2J_, all assumptions 
Tl) - (VIII) are needed. 



Some of the properties specified in Hypothesis |4.1[ 4.3, 4.4 and 4.5 have already got 
justified in Subsection |4.1 Moreover, it follows readily that P is affine. It remains to 
prove that P satisfies Hypothesis 4.3 c) while R satisfies Hypothesis 4.1 b). 



Lemma 4.6 If v' ^ v in L°°(0,T; 



K in K and = P(f*) then n = 'P{v). 



The proof of the above lemma follows along the lines of part of the proof of [TT| theo- 
rem 4.1]. 

□ 

Lemma 4.7 If C ^ C G Ms, then the sequence /* = R(^') has a subsequence 



convergent in sense of (4-3) to f = R(0- 



Proof. As in the proof of [TT| theorem 4.1], we are going to show that the sequence 
{f'^)°^i has a subsequence such that for every k = 1, . . . ,n for every 6k > and 

t G (0,T) the inequality | + h) — fk{t) | < 6k holds for s large enough (s > So{t,6k)) 
and h small enough {0 < \h \ < hQ{6k) with hQ{6k) G (0, 6k]). 

Consider the sequence for some k G {l,...,n}, where /* are defined as in 

the statement of the lemma. By the definition of the operator R and by the triangle 
inequality, 

\flit + h)-fkit)\ = \uie)ixlt + h)-uiOixl,t)\ 

< \u{e){xit+h)-u{e){xit)\ + \u{e){xit)-u{e){x,i)\ 

+ \u{C){x,i)-u{0{x,i)\ + |u(0(5,t)-u(e)(4,t)| (4.4) 

holds for arbitrary [x, t) G Qt- We will individually estimate each of the above four right 
hand side terms. 



By Theorem 5.3, the first of these terms can be estimated by: 



u{e){xi,t+h)-n{e){xi,t)\ < c,\\\{xi,t+h)-{xi,t)\\r 

= C6 111(0, /i)| I r = C6(V4)" <5fe/4 (4.5) 



is defined in the statement of Theorem 



5.3). 



for h < 4(4/(4c6))^/" (where 
Theorem 



5.3| allows also to infer analogous estimates for two of the remaining three 
terms in (4.4): 



u(r)(x*,t)-u(e)(x,t)| < C6 |||(x* -x,t-t)|||" < 6k/4: (4.6) 
u(0(x,t)-u(0(4,t)| < ce\\\ii-xl,i-t)\\r < 6k/4: (4.7) 



for {x,t) close enough to (a;^,t), i.e. for 
{x,i) e B := {ix,t) : |||( 



X 



Xki ^ V 



<(4/(4c6))^/"} . 
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Now we are left to estimate the remaining term in (4.4), namely | u(^*)(5, i) — u(^)(x, i) |. 
To this purpose, we will show that the sequence (u{^^){x, i))^^ has a subsequence, denote 
it (u(^*'')(a;, such that, for some {x,i) E B, 



u(^*^)(x,t) u(^)(x,t) ass-)-oo. 
By Theorem 5.2 and assumption ( |VIb ) : 



(4.8) 



u(e)-u(OL.(Q^) < C5 II <?(.,. ;e(. ))-(?(.,. ;e(.)) 

m 



Ili(Qt) 



I C[0,T] 



< cscsj^ll^ 
i=i 

the latter right-hand side is convergent to by assumptions of the lemma, hence there 
exists a subsequence (u(^*=))^]^ convergent to u(^) a.e. on Qt- In particular, this sub- 
sequence is convergent a.e. on B and hence we can choose {x, t) E B such that the 
convergence (4.8) holds, i.e. 



u{e')ix,i)-u{0{x,i)\ < 5k/ 4: 



(4.9) 



for s > So{Sk) and, besides, (4.6), (4.7) are satisfied . T hu s, co n side ri ng th e subsequence 
(i^(C''))^i aiid the point (x,t) as above and taking (4.4), (4.5), (4.6), (4.7) and (4.9) into 
account, we obtain 

\rkit + h)-Mt)\ < 4 

for s > So{Sk) and h < A ((5jt/(4 Ce))^/" what concludes the proof. 

□ 

To sum up, Ms is a nonempty, convex and compact subset of the Banach space K 



of admissible controls, the operators P, Q and R fulfill Hypothesis 4.1, 4.3, 4.4 and 4.5 



thus due to Theorem 4.2 we can deduce that the system (1.2) has at least one solution 
in sense of Definition 2.3 (or, equivalently. Definition 4.1) given by (u(/t), k), where k is a 
fixed point of P o Q o R. 



5 Appendix 

Below, technical definitions and results utilized in other sections of present paper will be 
provided, in part originating from other publications (P^, [7], [T5], [TS]). Still some of the 
formulations are adapted to the context of our paper, thus require more comments. 

Definition 5.1 For two arbitrary Banach spaces X and Y , a multivalued function F : 
X — 2^ is said to he upper semicontinuous if and only if, for every x E X and for every 
open set O C Y such that F{x) C O, there exists a neighborhood U{x) of x such that 
F{U{x)) C O. 

Remark. Note that for singlevalued functions the above definition reduces to the 
definition of continuity, not just semicontinuity. Thus, the introduced property is stronger 
than the standard semicontinuity. 
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Definition 5.2 For two arbitrary Banach spaces X and Y , a multivalued operator T : 
X — 7- 2^ is said to he closed if and only if, whenever Xn ^ x in X , yn ^ y in Y , 
y-n G T{xn) then y G T{x). Equivalently, we can say that the graph of T is closed in 
X xY. 



Definition 5.3 For a Banach space X and a multivalued operator T 
element x & X is said to be a fixed point ofT if and only if x ^ T{x). 



X 



an 



Theorem 5.1 Suppose that assumptions I^IIIa ), (Illb), (Via), (VII) are fulfilled. 
Then, for arbitrary k, & K the weak solution u of problem ( 2.3) in sense of Definition 2.1 
exists, is unique and belongs to L°°{Qt)- 

Proof. The proof of the existence and uniqueness can be done with the methods similar 
to those in the proof of [13^ Theorem 8.4] (still, the function spaces and restrictions for / 
there considered are slightly different from the ones assumed in |13| what should be taken 



into account when showing the suitable estimates for the proof of Theorem 5.1 ). A formal 



idea of the boundedness proof consists in showing for a.e. x G that u{x^ t) is uniformly 
bounded in L°°{0,T). To prove the estimate for given a; G fi, we test the main equation 
of (2.3) with u multiplied by a smooth cut-off function ^ such that 



^e,h = 1 
^e,h = 



h' 



on B{x,e) x (0,T) , 

on {B{x,e + h) x (0,^))^ 

on {B{x, e + h)\ B{x, e)) x (0, T) . 



In the obtained integral identity we pass with h —> and then, using the Lebesgue 
differentiation theorem, with e — )■ 0. After that, the required estimate can be concluded 
by the Gronwall inequality. 

□ 



Theorem 5.2 Suppose that the assumptions of Theorem 5.1 are satisfied and Ui and U2 



are solutions of the system (2.3) corresponding to controls Ki and K2, respectively. Then 



there exists a positive constant c^, independent on the control and initial conditions, such 
that: 



II «1 - «2 WmQT) + II ^1 - "2 |li2(Q^) < 05 \\ g{ . , . ; Ki{ .)) - g{ . , . ; K2{ .)) 



Proof. After transforming the integral identity in Definition 2.1 to the form with all 
derivatives put on a test function, the proof can be performed with arguments analogous 
to those used in the proof of [121 Theorem 4.3]. 

□ 



Theorem 5.3 Suppose that assumptions of Theorem 5.1 and assumptions (IIIc), ( Vb) 
are satisfied. Denote also 



ixA) 





^2 


t 




• 5 -^di 


4 



for every {x,t) G IR"'^"'^. Th 



en: 
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a) For arbitrary k, ^ K there exist constants a G (0, 1) and Ce > such that 



\u{Xi,ti) -U{x2,t2) I < C6|||(Xi,ti) - (X2,t2)|ir 

for all {xi,ti), {x2,t2) G intQ^ (note that we exclude the choice of {xi,ti),{x2,t2) G 
dQr)- 

b) The constant a is not dependent on k, and the constant Cq depends on\\ g{. , .;/«(.)) ||^g(o T-LP(n)) 
continuously. 



c) If, additionally, assumptions (Va) and (Vic) are fulfilled and the choice of k is re- 
stricted to the set Ms, then the constant Cq can be chosen such that is not dependent 
on K, i.e. Cq = Cq{k,) is a constant function on Ms. 



Proof. A proof of [a]) can be found in [T], Theorem 4], however we still give here short 
consistency remarks. 

The definition of weak solutions to a parabolic equation considered in [Ij (whose 



particular case is (2.3)) differs from Definition 2.1 (the former assumes that the time 
derivative is put on a test function and that the test functions are of class C^{Qt) H 



{(f){T, . ) = 0}), hence first we need to transform the integral identity in Definition 2.1| to 
the form as in \T\. 

Besides, in [Ij the nonlinear term is assumed to satisfy /(s) < Cy | s | with respect to 



the values of solution, differing from our case (since assumption (III) allows e.g. f{s 



s + s). But due to Theorem 5.1 for a given initial condition Uq and a control k any 



solution u of (2.3) is bounded. It can be shown that a bounded solution u of (2.3) is also 



a solution of an analogous system with the nonlinear / term replaced by 

+C7 for s s.t. s > C7 
fis) = { f{s) for s s.t. s G [-€7,07] , 
— Cy for s s.t. S < — C7 

where Cy = ess sup^^^^-^^q^ \u{x,t) \ . / as above belongs to the class of nonlinear terms 
admissible in [T]. 

After having overcome the above difficulties, flT, Theorem 4] can be applied. 

The assertion |b]) can be deduced by a technically extensive analysis of the proof of 
[H Theorem 4]. The latter proof consists in using a Harnack-type inequality, shown by 
suitable cut-off techniques. 

The assertionjcj) is a consequence of |b|. Since Cq depends oii\\ g{. , .;/?(.)) r-Lp(n)) 
continuously and me latter value is bounded with respect to k for k G Ms, the values of 
Cq also are bounded and thus it suffices to choose a constant cg = sup^g^j^ Cq in [a]). 

□ 



Theorem 5.4 Let W be a multivalued function satisfying assumptions (IVa), (IVb) and 
(IVc), f G C((0, T); M") (where C((0, T); M") is considered with standard supremum 
norm) andW{t) := W{t,f{t)) on [0,T]. Then there exists a measurable selection, i.e. a 
measurable function V : [0,T] — > R, such that V(t) G W(t) for a.e. t G [0,T]. 
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□ 



The proof of Theorem 5.4 follows in the same manner as that of ^ Lemma 3.4]. 



Theorem 5.5 IfW, f andW are as in Theorem 5.4. then there exists a solution to (2.4) 



in the sense of Definition 2^ 



Proof. This is a direct corollary of Theorem 5.4 Since W obeys assumption (IVc), W 
is bounded and hence any selection V(t) G W(t) (for a.e. t G [0,T]) is bounded as well. 
By Theorem 5.4, there exists a measurable selection. A selection V G L°°{0,T) which is 



bounded and measurable belongs to L°°(0, T). The existence of a solution to (2.5) in the 
Caratheodory sense is known (cf. |H Chap. 3, Sec. 3]). 



□ 



6 Concluding remarks 

The methods introduced in this paper can be applied to a broad class of reaction-diffusion 
models, which share certain key properties with our model, including a whole range of 
AUen-Cahn models. These properties are: existence of solutions, their uniqueness and 
regularity and the system stability with respect to perturbations of a control term, as in 
corresponding theorems in Section |5| 



Once the existence question is answered for (1.2), one can pose a problem of optimal 
localization of control devices and measurement points. This is a direction of our present 
research. In the case of Lipschitz continuous Wj in (1.2) we expect to obtain a result 



which says that for a given reference state of the system and a given number of the 
control devices and measurement point it is possible to localize those devices and points 
in a way which is optimal with respect to the task of preserving the evolution of the 
system as close as possible to the reference state. Considering the latter restriction on Wj 
also allows us to formulate the uniqueness result complementing the existence theorem 
presented in this paper. 

An optimal implementation of the control algorithms will be included in Grzegorz 
Dudziuk's Ph.D. thesis (under preparation). 
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